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1. Introduction. When a family of classical paths focuses, the envelope of the trajec- 
tories forms a focal region called caustics. The non-existence of paths or a coalescence of 
trajectories are interesting classical phenomena pertaining to caustics, which have been 
studied extensively in various fields of physics (see, e.g., [1]). Their consequences on quan- 
tum physics are, however, somewhat involved. This is mainly due to the fact that the 
powerful and commonly used semiclassical approximation scheme [2, 3] breaks down when 
caustics occur, since in this approach the transition amplitude is obtained by summing up 
fluctuations around (a finite number of) classical paths which are assumed to exist. So far 
there have appeared numerous works on the problem of quantum caustics. On one hand, a 
generalized prescription for the semiclassical approximation in the path-integral framework 
has been devised [4, 3], and on the other, a detailed quantum analysis for caustics in the 
presence of higher order terms (in the semiclassical expansion) — which smear out and so 
defuse the singularity at caustics — has been pursued [4, 5] (see also [3, 6, 7, 8]). 

In this letter, we wish to provide a simple method to derive a closed form of the 
transition kernel on caustics in the semiclassical approximation. We also mention some 
novel quantum aspects observed on and around caustics, namely, the concentration of 
amplitude on the focal point, observed in the Gaussian slit (gedanken-)experiment [9]. For 
simplicity we restrict ourselves to one dimensional systems without higher order terms 
focusing on systems governed by the quadratic Lagrangian of the form 3 

L[x] = l -x 2 - l -\(t)x 2 -n(t)x, (1.1) 

for which the semiclassical approximation is known to be exact. The transition amplitude 
between the two given endpoints a, b e IR during the interval [0, T] reads [10, 11] 

e ^I[x cl \-ifm{\) ^ ^ 12 ) 

where m(A) is the Morse index of the harmonic potential A and I[x c \\ = J dtL[x c \\ is the 
action evaluated for the classical path x c \(t) connecting the two endpoints a, b. Clearly, 
when the classical path ceases to exist, the prefactor \d 2 I[x c \]/dadb\ becomes singular and 
so does the transition kernel (1.2). We shall see that the kernel at the singularity, i.e., 
on caustics, can nevertheless be obtained in the modified prescription [3] by making use 
of the unitarity relation of the kernel, which improves the previous treatment [5] based 

Note that any Lagrangian at most quadratic in position and velocity can be brought into this form 
by partial integration. We use the dot to denote time derivative x = dx/dt, and put the mass m of the 
particle unity. 
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on an abstract setting. Further, our analysis of the Gaussian slit experiment reveals that 
the quantum effect suppresses the susceptibility of the focal concentration against initial 
momentum fluctuations. 



2. Classical caustics. Let us first recapitulate the classical caustics (for a fuller account, 
see [12]). The equation of motion derived from (1.1) and satisfied by the classical solution 
x c \ takes the form, 



A\ x c \(t) = p(t) , where A x := - 
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(2.1) 



If we let x c \(p, t) be the solutions characterized by the initial momenta x(p, 0) = p, then the 
Jacobi field J(p,t) := dx c \(p,t)/dp is a solution of the homogeneous part of the equation 
of motion (2.1), 

A x J(p,t) = 0, (2.2) 

satisfying J(p, 0) = 0. Since the Jacobi field so defined gives a measure for the variation at 
later times under initial momentum fluctuations, we see that caustics (at t = T) occur when 
its final value vanishes J(p, T) = 0. Note that, on account of the relation d 2 I[x c \\/dadb — 
— 1/J(p, T), the classical caustics implies quantum caustics as well, as seen in (1.2). Let 
u, v be two linearly independent solutions of the homogeneous equation (2.2) subject to 
the initial conditions u(0) = 0, -u(0) = 1 and v(0) = 1, i)(0) = 0. (In terms of the Jacobi 
field, we have u(t) = J(p, t)/ J(p, 0).) Then the general solution of the full equation (2.1) 
is given by 

x(t) = av(t)+/3u(t) + s(t) , (2.3) 

where s is a special solution of (2.1) satisfying s(0) = 0. The constants a, (3 are determined 
from the initial position x(0) and velocity i(0), respectively. 

When considering the quantum transition amplitude in semiclassical expansion, it is 
important to know whether to a given Dirichlet boundary condition 

x(0) = a, x(T) = b (2.4) 

there exists a classical trajectory. Suppose first there is no caustic at t = T, that is, the 
Jacobi field J(p,t) (and thus also u(t)) does not vanish at t = T. Then, by choosing 
a = a and (3 = (b — av(T) — s(T))/u(T) in the general solution (2.3), one obtains the 
desired trajectory which is clearly unique. However, when caustics occur at t = T (i.e., 
J(p 7 T) =0), then the final position x(T) turns out to be independent of the initial velocity 
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x(0), and while a = a, (3 is left arbitrary since u(T) = 0. Thus whatever the initial velocity 
may be, all paths starting from x(0) = a ends at the final point 



x(T) = v(T) ■ a + s(T) . (2.5) 

If the given value b does not match this unique value, then there is no classical trajectroy 
satisfying the Dirichlet boundary condition (2.4). 

This phenomenon can be considered as a special case of caustics in geometric optics 
(see e.g., [3]), and when it occurs we shall call the potential A critical, otherwise non- 
critical. For those critical potentials, it is convenient to introduce the constant 

*(A) := ^ , (2.6) 

where here v can be any solution of the homogeneous part (2.2) with v(0) ^ 0. It is easy 
to show that k{\) is indeed independent of v, and in view of (2.5), it gives the stretching 
factor during the period [0, T]. The final point x(T) in (2.5) is called the conjugate point 
to a or focal point. 

If A is critical, then besides the stretching factor k(X), a further characteristics is given 
by the Morse index, which we now briefly outline. Let I[x] be the action functional for 
paths x with fixed boundary condition (2.4). The second variation of the action 5 2 I around 
a given classical path x c \(t) (which in case of critical potential A is assumed to exist) is a 
symmetric bilinear functional of the perturbations around x c \, and its kernel is given by 
the operator A\ in (2.1). From the consideration of Sturm-Liouville problem the operator 
A\ is known to possess a complete set of orthonormal eigenfunctions u n (t) [13, 14]: 



A\u n (t) = 

with 
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u n (t) = E n u n (t) , (2.7) 



«n(0) = u n (T) = 0; / dtu n (t) u m (t) = 5 nm . (2.8) 

The index of the bilinear functional (Hessian) 5 2 I is given by the number of negative 
eigenvalues of eigenvalue problem (2.7) for non-degenerate A\ for which no zero mode 
exists. It characterizes the type of the saddle point of the action functional I[x] at the 
extremum path x c \. (In this letter, we assume the eigenvalues to be bounded from below, 
which is for example the case when A is bounded from above.) The index is known to be 
equal to the Morse index which is the number of zeros of the Jacobi field on the half-open 
interval (0, T] [15] and denoted by m(A). Note that the bilinear functional A\ becomes 
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degenerate if there arises a zero mode solution of (2.7), a u m with E m = for some m, i.e., 
when caustics occur. In this case the Morse index gives the number of negative eigenvalues 
plus one, and accordingly the index of A\ may be extended even to this degenerate case 
by saying that it is given by the number of modes with E n < 0. 

3. Quantum caustics. We now study the quantum dynamics of the system (1.1). Let I[x] 
be the action evaluated along a path fulfilling the boundary condition (2.4). Corresponding 
to this boundary condition (2.4), the transition amplitude in the path-integral is given by 

rx(T)=b 

K(b,T;a,0)= / Vxe* 1 ^. (3.1) 

Jx(0)=a 

In carrying out the path-integration, we need to take into account the fact that, when 
caustics occur, there may not exist a classical solution that respects the given boundary 
condition. The general case, including the caustic one, may be handled by the following 
procedure [3, 4]. 

Let x c \ be a classical solution of the equation of motion (2.1) with the initial value 
^ci(O) = a, and let c := x c \(T) be the endpoint of the solution. Note that in general c may 
not be equal to 6, but if there is a solution having b as its endpoint, then we shall omit 
the bar and denote this special solution by x c \. Such a solution exists if A is non-critical, 
whereas if A is critical the endpoint is determined uniquely c = k(X) a + s(T), see (2.5) 
and (2.6). Now, let us introduce a fixed, smooth function p{t) satisfying p(0) = and 
p(T) = b — c. Then any path obeying the boundary condition (2.4) may be decomposed as 

x(t) = x c i(t) + p(t) + rj(t) , (3.2) 

with rj{t) being the function representing the fluctuations vanishing at both of the ends 
i](0) = f](T) = 0. With the help of the orthonormal eigenf unctions in (2.7) and (2.8) it 
may be expanded as i](t) = J2 n anUn (^)- U s i n g this and the decomposition (3.2) we find 
by partial integrations 

I[xd + P + V}= I[x c i] + I[P\ U=o + *ci(T)p(T) 

i r T 

+ -J2 E nal + p{T)J2 a nUn{ T ) + J2 EnCln / & p(t)u n (t) 
n n n ^° 

In order to perform the path-integral (3.1) we change the integral variables Vx = Vrj = 
A/"n n ^ a n where Af is a Jacobian factor. Further, we replace E n with E n + ie, e > 
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(3.3) 



infinitesimal, such that the integrations over a n become all Gaussian even if an eigenvalue 
E m should vanish. With these preparations one obtains 



K(b, T; a, 0) = A/" 



Y[(E n + ie) 



(3.4) 



where 

$(&, a; A) := I[x d ] + I[p] | M =o + S d (T)p(T) 



(3.5) 



The kernel formula (3.4) is valid for any potential A. Now for A non-critical, each E n ^ 0, 
and in the limit e — > the kernel (3.4) reduces to 



K(b,T; a,0) = M 



11 E, 



(3.6) 



where p(t) has been set identically zero in (3.5) owing to the existence of the solution x c \. 
The obtained kernel (3.6) can be shown [10, 11] to be equivalent to (1.2). 

If, on the other hand, X(t) is critical, then E m = for some m, and the integration over 
a m yields in the limit e — > the delta- function S(p(T) u m (T)) on account of the identity, 
lim e ^o ~m — e ~ x ^ 2e = S(x) . Hence, for critical X(t) the kernel formula (3.4) reduces to 



K(b,T;a,0) 



2tt 




\{E n 



8(p(T)u m (T))e^^ 



(3.7) 



In (3.7) we have set p(t) = in the phase part of the kernel (and hence we put I[x c \] = 
/[x c i]) , which is allowed due to the delta- function and the fact u m {T) ^ 0. (This follows 
since both u m {T) = and u m (T) = would imply u m = identically.) We note that 



p(T) = b - k(X)a - s(T) 



(3.8) 



and, therefore, the kernel (3.7) clearly shows that classically forbidden transitions to non- 
conjugate points, b 7^ k(X)a + s(T) from a, remains forbidden even quantum mechanically. 

The remarkable fact in the case of caustics is that one can express the transition kernel 
(3.7) in terms of the stretching factor k(X) and the Morse index m(X). To see this, let us 
write the kernel (3.7) in the polar form, 



K(b,T;a,0) = R(T)S(p(T))e i ^ T \ 



(3.9) 



6 



Noting that the transition kernel K(b, T; a, 0) can be expressed as (6| U{T, 0) \a) in terms 
of time evolution operator U (T, 0) , the unitarity relation, 

J dbK*(b,T;c,0)K(b,T;a,0) = J db (c\U\T,0) \b){b\U(T,0) \a) = 6 (a - c), (3.10) 

determines the modulus R(T) to be R(T) = y/\k(X)\. On the other hand, since the 
Morse index gives the number of negative modes plus one, we have [Iln^m -^] _ ^ = 
[rin^ml^nr^e-^^)- 1 ). Thus the phase part is given by 6(T) = J/[x c i]-f (m(A)-l). 
By combining these we find that the transition kernel (3.7) takes the closed form in terms 
of the stretching factor and Morse index, 

K(b, T; a, 0) = ^/\kjxj\ S(b - k(X)a - s(T)) e * /[Xcl] -^ m(A) , (3.11) 

up to an overall constant of unit modulus. This final form agrees with the expression 
obtained earlier based on a different construction of the path-integral [5]. 

As application, take, for example, the forced harmonic oscillator given by X(t) = uj 2 
and n(t) = —f(t). A special solution satisfying s(0) = is then 

1 f* 

s (t) = - / dt' sinw(t - t') f(t'). (3.12) 
u Jo 

For the present case, caustics occur at u = ^ with n = 1,2, . . ., where we have k(X) = 
(— l) n and m(A) = n, respectively. The classical trajectory starting from x c \(0) = a and 
ending up at the point (2.5) possesses the action, 

I[x cl ] = a f dt cos utf(t) - - f dt [ dt' cos ut smut' f(t) f(t') . (3.13) 
Jo v J J 

Hence, in particular, for constant / the kernel (3.11) reduces to 

K(b, T; a, 0) = s(b - (-l) n a - {1 - (-l) n }//w 2 ) e if 2 T/2^^nn/2_ (314) 

This formula coincides with the previous results [16, 17, 18] obtained by other indirect 
means 4 for the (forced) harmonic oscillator. 

4. Gaussian slit experiment. For A non-critical but close to a critical A, one expects 
that there will be a concentration in the transition amplitude around the focal point 

4 The result in [18] obtained for f(t) ^ is marred by an error in the phase. 
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b = k(X)a + s(T) conjugate to a given initial point a. We analyze how the concentration 
takes place quantum mechanically, based on the Gaussian slit (gedanken) experiment [9]. 

Let there be an apparatus which emits a particle from the origin x = at time t = —t. 
To get a Gaussian distribution at t = 0, we place a 'Gaussian slit' centered at x = a with 
effective width (variance) o"o, which is given by the Gaussian transmission factor 

G(x - a; <T ) = N exp j- ^"^ j . (4.1) 

The wave function ip(x, 0) of the particle at t = is then furnished by the product of 
the free particle kernel K (x, 0; 0, — r) = ^J1/(2-kiTit) e lx2 / 2hT and the transmission factor. 
From the normalization condition f dx \ip(x, 0) | 2 = 1 we determine the constant iV and 
obtain 



1 J (x — a) 2 ix 2 



V>(x, 0) = G(x - a; (to) K (x, 0; 0, -r) = - exp { - V - 2 J + — \ , (4.2) 



which has the average momentum p = J dxip*(x 1 0) [—ih-^)ip(x 1 0) = ajr. Using the 
kernel (1.2), the wave function at t = T is then given by 

V> (y, T) = J dx K(y, T; x, 0) ^(x, 0). (4.3) 

To simplify the ongoing discussion let us assume first /j, = in (1.1). Then, the classical 
action 7[x c i](y, T; x, 0) is a quadratic polynomial homogeneous in x and y, because we have 
I[x c \](cy, T; cx, 0) = /[cx c i] = c 2 /[x c i] = c 2 /[x c i](y, T; x, 0) for any constant c. Thus, with 
A, S and C some functions of T, we may write 

7[x c i] (y, T; x, 0) = Ax 2 + Bxy + Cy 2 . (4.4) 

Then it is straightforward to see that the wave function (4.3) at time t = T has the form 

if>(y,T) = r exp I - (y ~ ^ + i (phase) ) ( 1.5) 

VKy ' (2tt(7 2 (T))i l \ 4a 2 (T) ^ ^ 



with 



and 



x cl (T) = -^(2aA + p), (4.6) 
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The second term in (4.7) represents the quantum effect whereas the first term is just the 
classical variance, since the variation in the initial position Ax c \(0) = o"o around the center 
^ci(O) = a for the classical path results in the final variation Ax c \(T) = a \x c \(T)\/a = 
cro\k(X)\. (Note that the initial momentum p is linearly dependent on the initial position 
where the path goes through.) Thus, as one expects, the quantum effect always enhances 
the spread of the Gaussian distribution. 

Viewed as a function of the initial variance o"o, the final variance cr(T) attains its 
minimum, i.e., the highest concentration of intensity 



hx c \{T) 



aB 



at Co = 



ha 



2Bx cl (T) 



(4.8) 



which is precisely the point where the quantum effect matches the classical variance. It 
then follows from (4.6) and (4.8) that, ideally, the infinite concentration cx m i n (T) = takes 
place at B = ±oo or at p = —2aA. The former case corresponds to the expected caustics, 
which occur when the variance (4.7) becomes purely classical cr(T) = ao\k(X)\ for which 
the infinite concentration is obtained by letting o~o — > 0. By contrast, in the latter case 
the variance becomes purely quantum mechanical, as all paths passing through the slit 
coalesce toward the origin x = at t = T. In fact, this case is again caustics occurring 
during the combined period [— r, T] under the potential X(t) which vanishes for t < 0, 
where the infinite concentration at the quantum level is achieved if we let <jq — > oo. It is 
obvious, however, that both of these concentrations in intensity are unstable, since a small 
fluctuation in the parameters B or p will generally bring the variance to a large value. 
Nevertheless, one could achieve a sharp peak in intensity by adjusting the initial variance 
along with the parameters according to (4.8). 

An important quantity for characterizing the concentration is the susceptibility of the 
variance against initial momentum fluctuations. Exposing the momentum dependence of 
the variance explicit cr(T) = a(p,T), we use the following normalized quantity for the 
susceptibility, 

The susceptibility of (4.7) reads 

*" = '* T »f(^m) '> ■ <410) 

where J{p,T) = — [d 2 I[x c \]/ dydx) 1 = —1/B is the Jacobi field for the classical action 
(4.4). In the classical limit % — > 0, the susceptibility reduces to the (absolute value of the) 
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Jacobi field. In the two cases mentioned above in which an infinite concentration can in 
principle be possible, the susceptibility vanishes S(p, T) = and hence the Gaussian wave 
packet becomes free from momentum fluctuations. What is interesting in (4.10) is that 
the quantum effect suppresses the susceptibility of the variation against initial momentum 
fluctuations. 

The above features of the Gaussian wave packet persist even for /i^O, because then 
the classical action I[x c \](y,T;x,0) acquires only linear and constant terms in addition to 
the quadratic terms in (4.4) as can be explicitly seen from the form of the general classical 
solution derived after eq. (2.4). Accordingly, the time evolution by the integral (4.3) is 
essentially unchanged. Semiclassically, this may also be the case for more general systems, 
not only for those with quadratic Lagrangians we considered, in view of the earlier study 
[19] which suggests that these features are a norm for a generic system in the limit % — > 0. 

As an example of the above discussion, consider again the harmonic oscillator A = uj 2 
(but with /j, = 0). ^From the familiar classical action [9] one reads off that A = C = 
uj/2taiYUjT and B = —uj/ sinwT. With the center of the distribution x c \(T) = a cos lvT + 
ip/ui) smivT, one finds the variance, 

^)-{(co S . r+ £ S in.r)%(^) 2 }\ (411) 

possessing its minimum cr m i n (T) = \hx c \(T) smu)T/auj\^ at do = | ha sin u>T/2u>x c i(T)\^ . 
One thus sees that an infinite concentration of the amplitude takes place either at u = ^ 
or p = —au>/ tancuT. 

5. Conclusion. When caustics occur in one dimension, only those paths connecting the 
conjugate points are allowed classically, and the intrinsic features of caustics are character- 
ized by the stretching factor k(X) and the Morse index m(A). In quantum mechanics the 
transition is still allowed only between conjugate points, and for systems with quadratic 
Lagrangians we derived the path-integral kernel for the amplitude in a closed form, ex- 
pressed solely in terms of the stretching factor, the Morse index, and the action of (any 
of the) solution paths. As for the situation near caustics, the Gaussian slit experiment 
showed that, although the variance itself is enhanced at the quantum level, the suscepti- 
bility of the variance of the packet against initial momentum fluctuations is suppressed. 
This was derived by considering a quantum analogue of the Jacobi field, which reduces to 
the ordinary Jacobi field in the classical limit. Possible extensions of this work on quantum 
caustics, e.g., effects of higher order terms in the Lagrangian on the quantum susceptibility 
or related problems in higher dimensions, are currently under investigation [20]. 
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